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A continuum of structure and stellar content from Virgo cluster early- 
type dwarfs to giants? 
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Based on the wealth of multiwavelength imaging data from the SDSS, we investigate whether dwarf and giant early-type 
galaxies in the Virgo cluster follow a continuum in their structural parameters and their stellar population characteristics. 
More specifically we study the relation between size and brightness for the galaxies and their color magnitude relation. In 
both cases, we find noticeable deviations from a simple joint behavior of dwarfs and giants. We discuss these findings in 
the light of the different formation mechanisms commonly assumed for dwarf and giant early types, thereby taking into 
account the existence of several distinct early-type dwarf subclasses. By comparing our results to a semianalytic model 
of galaxy formation, we argue that the analyze d relations mi^ht be reprod uced by processes that form dwarfs and giants 
altogether. The work presented here is based on ljanz & Liskei 



ignt be reprodi 
ej|2008l , l2009l 



1 Introduction 

Early-type dwarf (dE) galaxies are commonly expected to 
play a key role in understanding galaxy cluster evolution. 
Their importance is given by their abundances - they out- 
number all other galaxy types in dense cluster environments 
by far - and the fact that they provide not too massive, not 
too dense test particles to probe processes that let the cluster 
environment alter the appearance of galaxies. At the same 
time, dEs are predicted to form in models of a ACDM uni- 
verse as the descendants of building blocks in hierarchical 
structure formation and to be in that sense close relatives 
to their giant counterparts, sharing a cosmological origin. A 
better understanding of dEs is therefore not only linked to 
our knowledge of formation and evolution of galaxy clusters 
but also of structure formation itself. 

Once believed to be systems of simple appearance and 
well-defined properties, dEs were recently shown to exhibit 
a puzzling variety among their structures and stellar popula- 
tions (see e.g. T. Lisker, this issue). This diversity opens the 
door widely for different formation scenarios. And indeed 
there are different suggestions, for example the transforma- 
tion of other galaxy types by the cluster environment via 
ram pressure stripping or harassment, which are partly able 
to explain the appearance of dEs and also reproduce with 
some successes fundamental scaling relations of early-type 
galaxies. But still today, it remains an open question to what 
extent these different processes play a role and whether some 
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of the early-type dwarf galaxies share the same origin and 
formation mechanisms with their more massive relatives. 

The above mentioned scaling relations have ever been 
an important tool not only to study galaxy properties but 
also to link those properties to their formation and evolu- 
tion, and thus to answer the question. Very well studied 

examples are the relations between surface brightness and 

I II — — I 

size ("Kormendy relation". IKormendvf 1985). betwee n sur- 
face b rightness and luminosity (e.g. IBinggeri & Cameronl 
Il99lh . In combi nation with velocity disp ersion, the Faber- 
Jackson relation dFaber & Jacksonlll 976l) and the e xt ension 
the Fundamental Plane dDressler et al.1 11987 , 



to 



Diorgovski & Davisl 1987b became famous. Every time these 



relations were analyzed for dwarfs and giants in conjunc- 
tion, it was discussed, whether or not they show a common 
behavior and what causes it. Any dwarf formation scenario 
has to reproduce the observationally found relations. Addi- 
tionally to these morphological and kinematical relations, 
the color magnitude relation (CMR), connecting the global 
parameter total brightness of a g alaxy to its stellar popula- 
tion, w a s exten s ively s t udied. e.g.[B^uml(ll959l):lBpwer et al.l 
dl992b: iFaberl dl973h: ISandaee & Visvanathanl dl978albF: 
Visvanathan & Sandagd d 1977b . The CMR is typically ex- 



plained bj^anjiicre^eofmean stellar metallicity (and age; 
see e.g. Gavazzi et al . 2002) with increasing galaxy mass as 
the dominant effect. The common underlying idea is that 
more massive galaxies have deeper potential wells, which 
can retain metal-enriched stellar ejecta more effectively and 
subse quently recycle the enriched ga s into new 
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stars (Chang et al. . 2006 : Ferreras et al. . 19991 Gallazzi et al 
2006; Kodam a & Arimotoll 1997b " Also here it was explored 



'*WILEY 

jinterScience* 



© 2009 WILEY-VCH Verkg GmbH & Co. KGaA, Weinheim 



2 



J. Janz & T. Lisker: Dwarf and giant early types 



whether and, if so, how much these processes shape the 
CMR of giants and dwarfs in a similar way. 

We made use of a very homogenous data set of the early 
types in the Virgo cluster to investigate these questions via 
the scaling relation of size and brightness (which is a 
relative to the aforementioned morpho logical scaling rela- 
tions) and the color magnitude relation dJanz & Liskerll2008 
2009h . 



and j = 0, . . . , 16. We omitted intensities at radii r < 2" in 
order to avoid seeing effects. 

Colors were measured within the elliptical r-band half- 
light aperture for each filter. Errors were estimated from 
the S/N and calibration uncertainties (which we estimate 
to have a relative effect of 0.01 mag in each band, which 
is smaller than the absolute v alues given by SDSS), as de- 
scribed in lLisker et al. (2008). 



2 Sample Selection and Imaging Data 



3 Sizes of early-type galaxies 



Our sample is based on the Virgo Cluster Catalog (VCC; 



Bing geli et al.l ll985). All early-type galaxies therein with 
a certain cluster member status and ms < 18.0 mag are 
taken into account, which is the same magnitude limit up to 
which the VCC was found to be complete. This translates 
into Mb < — 13.09 mag with our a dopted distanc e modulus 
of m-M=31.09 mag (d=16.5 Mpc. lMei et alll2007h . 

Uncertain classifications are treated as follows: galaxies 
listed as "SO:", "E/S0", "SO/Sa", and "SBO/SBa" are taken 
as SO, and one SO (VCC 1902) is excluded, since it shows 
clear spiral arm structure. For the dwarfs, we selected galax- 
ies classified as dE, dSO, and "dE:", whereas "dE/Im" as 
well as possible irregular s based on visual inspection are ex- 



rreg 

eluded dLisker et all 120071) . We exclude 37 galaxies for the 



following reasons: the Petrosian aperture (see below) could 
not be obtained, the objects were too strongly contaminated 
by the light of close neighbour objects, or the S/N in either 
the u or the z band was too low. Our working sample thus 
consists of 468 galaxies. 

The Sloan Digital Sky Survey (SPSS) Data Release Five 
(DR5) dAdelman-McCarthv et all 120071) covers all but six 
early-type dwarf galaxies of the VCC. Since the quality of 
sky level subtraction of the SDSS pipeline is insufficient, 
we us e sky-subtracted images as provided by iLisker et al.l 



(2007), based on a careful subtraction method. The images 
were flux-calibrated and corrected for galactic extinction 
dSchfegeletallll998l) . 

For each galaxy, we deter mined a "Petrosian semima- 
jor axis" a p dPetrosian , Il976l) . i.e. we use ellipses instead 
of c ircles in the cal culation of the Petrosian radius (see, 
Lotz et al 4). total flux in the r-band was mea- 



e-g 

sured within 2a p , yielding a value for the half-light semi- 



3.1 Introduction 

The scaling relation of size and brightness of early types 
was not as widely studied as its relatives, like the Kormendy 
relation and the relation between brightness and effective 
surfac e brightness. St udies of the sizes, for e xample , 
are: iBender et all dl992l) : iGuzman et all dl993l) : 



Kormendv l d 19771). and for th e Virgo Cluster in p articular by 



Binggeli & Cameronl d!991l) for dwarfs and by ICaon et al 



d 1993b for giants. But all of them share a similar history: 
Early studies until the 1990's came to the conclusion that 
giant and dwarf early-type galaxies show a distinct behavior 
in the scalings such as the relation between size and bright- 
ness. The dwarfs were seen to show less change of size with 
luminosity than the giants. This together with the other scal- 
ing relations was interpreted as evidence for a different ori- 
gin of dwarf and giant early-type galaxies. 

Towards the turn of the millenium, however, it became 
more widely realized that the light profile shapes of early 
types vary continuously with luminosity. Neither do dwarf 
galaxies simply follow exponential profiles, nor do all gi- 
ants exhibit de Vaucouleurs profiles. Instead, all earl y types 



are w ell described by the generalized Sersic profi le (Sersic, 
1963) with different Sersic indices n dFerrarese et all l2006: 
Young & Curriel 1 1 9941) . Several authors reasoned that the 



scaling relations naturally follow what is predicted by n 
changing linearly with magnitude, and that all these galax 



ies can indeed be o f t he sam e kind ( Jerien & Binggel 



Binggeli & Jerien 19981: Graham & Gu zman 



li 1 



1997 



2003 



Gavazzi et al.l2005l) . In ljanz & Liskerld2008l) we studied the 



size brightness relation of early types in Virgo and analyzed 
it in the light of a continuous variation of profile shapes. 



. This Petrosian aperture still 3 2 Results 



major axis (SMA), au,r,unaorr 

misses some flux, which is of par ticular relevance for the 
giant galaxies dTruiillo et alll200ll) . The brightness and the 



half-light SM A were correct ed for this missing flux accord- 
ing to iGraham et alj ( 12005b . Axial ratio and position an- 
gle were then determined through an isophotal fit at 2aya^ r . 
The effective radius is then given by r e g — a^j r \Jb/a 
with the axis ratio b/a. Additionally we fitted Sersic pro- 
files to the radial intensity profiles, keeping half light ra- 
dius fixed and using an implementation of the nonlinear 
least-squares Levenberg-Marquardt algorithm. For the fits 
we used the intensities at r/a^ r = 2 X with x = —2 + j/4 



In Fig. 1 (bottom panel) we present the size luminosity di- 
agram for our sample. At first glance the sequence from 
dwarf to giant early-type galaxies does not look very contin- 
uous: the giants follow a steep relation with a well-defined 
edge on the bright end of their distribution. The bunch of 
dwarfs apparently lie with a larger scatter around an effec- 
tive radius of r e ^ — 1 kpc, their sizes showing weak to no 
dependence on luminosity. 

A similar impressio n can also be obtained from other 
previous studies, e.g. Binggeli & Cameronl 1 1 99 ll Fig. lb; 
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Fig. 1 Bottom panel: Absolute magnitude in r versus log- 
arithm of half light radius. Filled squares - E, gray stars - SO, 
open squares - M32 candidates, open triangles - dE,N, open 
pentagons - dE,nN and gray open circles for dwarf galax- 
ies with probable disk-like structure (dE(di) or dE(bc)). The 
grey line is calculated with the linear fits in the top panels 
(see text). 



Bender et al.ll 19921 Table h lKormendv etal.ll2009l Fig. 37 



It is, however, not as clearly seen in the compilation of sizes 
of elliptical galaxies from sev eral different studies presented 
by iGraham & Worleyl (12008b (Fig. 10). In this more hetero- 
geneous data set, the relative number of small low-luminos- 
ity giants as well as that of large bright dwarfs appears to be 
somewhat smaller. 



Varying profile shapes ? 



Graham & Guz man (2003) suggested that the apparent di- 
chotomy between dwarfs and giants in scaling relations can 
be explained just by the fact that the profile shape of a galaxy 
scales with magnitude. They describe the light profiles with 
Sersic profiles and show the effect of a linear relation be- 
tween magnitude and logarithm of the Sersic index n on 
the other scaling relations. As a result, the dependence of 
effective radius on magnitude becomes stronger at higher 
luminositi es and the brightest galax ies are naturally larger 
(Fig. 1 1 in lGraham & Worlevll2008l) . 

For investigating whether our galaxies display the pre- 
dicted behavior, we use the Sersic indices n and central 
surface brightnesses /io to obtain linear fits to the fio /M r 
and n/M r relations, using a least squares fitting algorithm 
(Fig. 1, top panels). For those fits we exclude systems with 
a (probable) disk component , namely galaxies cl assified as 
SO, dEs with disk features ( Lisker et al., 2006a), and dEs 
with blue centers dLisker et all 12006b). This is to ensure 
that the light profiles can be well parametrized by Sersic 
profiles. Our fits (Fig. 1, top panels) together with equation 



(16) of IGraham & Worlevl (|2008) predict a non-linear se- 
quence in the r e j^/M r diagram. The predicted relation is 
shown together with the observed galaxies in the bottom 
panel of Fig. 1. With the visual guidance of the line, it ap- 
pears more likely that the data points follow one common 
continuous relation. And the gross trend in the diagram can 
indeed be explained by varying profile shapes. However, at 
luminosities brightwards of the transition between dwarfs 
and giants, a substantial amount of galaxies fall below the 
relation, while faintwards most of the dwarfs lie above it. 

As we showed in Janz & Lisker (2008) this finding holds 
also if all objects with signs of disk components are omitted 
in order to have a purer sample of dynamically hot systems 
not biased by systems with more complex kinematics. Fur- 
thermore, we analyzed the deviations statistically and found 
significance (Fig. 2 therein). We will discuss the implica- 
tions in Sect.[5]and Sect. [6] 

Our analysis showed two things. First, the distribution of 
data points does not resemble a quite large random scatter 
around the relation. And therefore the size luminosity re- 
lation can not be fully explained by varying profile shapes. 
Second, the abrupt change in the behavior of faint and bright 
galaxies is even emphasized through the above examination, 
and this break is a real discontinuity of the sequence from 
lowest to highest luminosities. 

4 Color magnitude relation 
4.1 Introduction 

From early on, it was discussed whether the universality of 
the CMR also holds over the whole range of galaxy masses, 
i.e. whether dwarf and giant early-type galaxies follow the 
same CMR. Studies of different clusters show consistency 
with one common CMR for dwarfs and giants, albeit with a 
significant increase in the scatter at low luminosities 



(Seek er et alj 1 19971 for Co ma IConselice et alj |2002| for 



Perseus, Karick et al. 2003 and iMieske et al. 
nax, Smith Castelli et al. 200£ 



2007 



for For- 



_.. . ._ _ . _ for Antlia, and Misgeld et al. 

2008) for Hydra I). More explicitely, ICaldwellld 1983b stated 



that there is a common linear relation. But his Fig. 3 might 
hint at a chan ge of slope from high to low luminosities, 
similar to what Ide Vaucouleursl dl96ll) suggested. Interest- 
ingly, visual examination of the diagrams presented by most 
of the above-mentioned studies indicates consistency also 
with a change of slope - yet linear relations we re fitted in 
most cases (see, however. iFerrarese et alj|2006t our colors 



are consistent with the ones in their Virgo c luster study in 
the ran ge of brightness common to both). In Ijanz & Liskei 
(2009) we revisited the question of the universality of the 
CMR for dwarfs and giants. 

4.2 Result 

We sh owed the CMRs for fo ur different representative col- 
ors in Jan z & Liskei (2009). Here we choose to show the 
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Fig. 2 Color Magnitude Relation in u — z. Filled cir- 
cles show our sample's galaxies. The gray line indicates the 
"running histogram" as found in successive magnitude bins 
with a width of 1 mag and steps of 0.25 mag, clipped one 
time at 3o\ We limit the drawing range for the line to the 
region with at least three galaxies in a bin. The white his- 
tograms show the distributions in bins of the same width, 
normalized to the square root of the number of galaxies in 
the bin, shown for every fourth step. The white errorbars 
indicate typical photometric errors at the respective bright- 
ness. In the panel right to the color magnitude diagram a 
measure for the intrinsic scatter is given by the RMS differ- 
ence S of the observed scatter and the photometric error (see 
text) in continuous bins. 



CMR in u — z, the color with the longest wavelength base- 
line available to SDSS. It looks remarkabl y similar to the 
CMRs in the other colors (Fig. 1 in I Janz & Liskerll2009l) . 

First of all, the impression one can get by examining just 
the black points in Fig. 2 is that there is not one common lin- 
ear relation from the faint to the bright galaxies. The overall 
shape appears more like "S" shaped. The brightest (M r < 
—21) galaxies have almost constant color, i.e. no correlation 
between color and brightness; the very brightest galaxies 
show a larger scatter. These were reported before to be mor- 
phologically different from the other galax ies in more de- 



tailed stu dies of the inner lig ht pro file s (e.g. Ferrarese et al 



20061: iKormendv et all 120091: lLaueretal.1 [2007 



Ravindranath et alj I200U iTruiillo et alj |2004 . For the 



maining galaxies several descriptions seem to be plausible, 
ranging from just an offset between two relations with sim- 
ilar slopes up to a curved relation. 

With the non-linear shape, it seems not very favorable 
to fit a straight line. This would not describe the data well, 
and there is no theoretical prediction what other function 
is expected. So at first, we want to make the overall shape 
more clearly visible, using continous, overlapping magni- 
tude bins, in which mean and scatter are calculated. In Fig. 2 



these derived relations are shown with grey lines. The first 
impression is confirmed: one common linear relation for 
dwarfs and giants cannot be seen. Moreover, the white his- 
tograms showing the galaxy distributions in the bins are 
clearly peaked towards the bright and the faint end, while 
they are rather flat at intermediate luminosity. 

The scatter about the relation is greatly influenced by 
increasing photometric errors towards faint brightness. In 
order to measure the intrinsic scatter, we calculate the RMS 
of the scatter around the mean in running bins (clipping one 
time at 3a) and subtract the RMS of the photometric errors 



S = RMS difference = 



s scat 



s err 



= jE(«-< c >i) a -E°?' 

V i i 

with color c and mean color (c), averaging over the galaxies 
in the respective bin. Here we exclude dEs with blue cores, 
since they are known to have different colors dlisker et al. , 
2008). This RMS difference should be zero if the scatter is 



only due to the measurement errors and larger if there is an 
intrinsic scatter. In Fig. 2 we show the CMR along with the 
RMS difference 6. 

Indeed, the RMS difference is enhanced for the dwarfs 
and peaks around M r « —18, indicating an intrinsically in- 
creased scatter, consistent with the increased intrinsic scat- 



ter for dwarfs found by See ker et alj (119971) in the Coma 
cluster. One can argue about the significance of the RMS 
difference increase for the brightest galaxies, since it is just 
a handful of them - nevertheless, this larger scatter might be 
related to the absence of a well-defined CMR at the bright- 
est magnitudes. 

5 Comparison to a Semi- Analytic Model 

Dark-halo merger trees of a high resolution TV-body simula- 
tion of ACDM structure formation were taken as input for a 
semi-analytic model (SAM) of the physical processes gov- 
erning galaxy formation and e volution in order to pro duce 



the Numerical Galaxy Catalog dNagashima et al.U2005l) . In 



particular, the dynamical response to starburst-induced gas 
removal after gas-rich mergers (also for cases intermedi- 
ate between a purely baryonic cloud and a baryonic cloud 
fully supported by s urrounding dark matter as in 



Yosh ii & Arimotoll 1987b were taken into account. This pro- 
cess plays a crucial role for the sizes of early-type dwarf 
galaxies, since the subsequent variation of the potential re- 
sults in an increase in size. 

We identify model galaxies as early-type galaxies if they 
are bulge dominated (bulge to total ratio > 0.6). To compare 
our data with the mo del we transformed SDSS g magnitudes 
into B according to lSmifh et al. ( 2002 ). using the galaxies' 
g — r color measured within ahi, r - B-V was obtained like- 
wise. In Fig. 3 one can see how the model galaxies com- 
pare with our observed Virgo early types. In the left panel 
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Fig. 3 Comparison to the semianalytic model by Nagashima et al. (2005). In grayscales (on a logarithmic scale) the 
distribution of model galaxies is shown. Model galaxies with fainter surface brightnesses than a limit of (/is) < 
25.5 mag/arcsec 2 are excluded. The observed Virgo early type are displayed with black open circles. Left panel: Ab- 
solute B magnitude versus logarithm of half light radius. Right panel: Color Magnitude Diagram. Colors and magnitudes 
were transformed to the filterbands of the model output (see text). 



the comparison is displayed for the size brightness relation. 
The model galaxies show a bimodality similar to what we 
observe, with low galaxy d ensity between the two regions. 
Note that Nag ashima et alJ assume de Vaucouleurs profiles 
to calculate projected half-light radii from half-mass radii. 
For exponential profiles, which would be more appropriate 
for dwarfs, the mode l galaxies would shift upwa rds in the 
diagram by 0.1 1 dex dNagashima & YoshiiLl2003l) . 



In the model, a starburst follows in those dwarfs that 
form by gas-rich major mergers and the dwarfs are enlarged 
by the dynamic response to the subsequent gas loss. This 
mechanism is not at work in gas-deficient mergers, and the 
resulting galaxies stay smaller. Note, tho ugh, that this ap- 



pears to be in contrast with the SAM of iKhochfar & Silk 



d2006l) . where gas-rich mergers lead to more compact early- 
type galaxies than gas-poor mergers. 

In the right panel of Fig. 3 the comparison is displayed 
for the CMR. The shapes of the distribution of Virgo galax- 
ies and the model CMR are indeed not well represented 
by linear relations. Both the CMR and the relation between 
size and brightness do not show a linear, nor one common 
behaviour of dwarfs and giants. Nevertheless, it needs to 
be emphasized that there is a qualitative similarity in the 
shapes of observed and model CMR, in the sense that both 
show a similar bend at intermediate luminosities. This is 
noteworthy, since in the framework of the SAM, both dwarf 
and giants form by the same physical processes, which gov- 
ern ACDM structure form ation, and thus both can be of cos- 
mological origin (see also lChilingarianll2009h . 



Beside the similarity in the overall shape, an offset is 
observed. This offset could partly be due to uncertainties 
of the adopted synthetic stellar population model. Further- 
more, the relative number of bright galaxies exceeds the 
observed one and the luminosity function is clearly differ- 
ent, which could possibly be explained with model input 
physics. 

6 Discussion 

We studied two scaling relations of the Virgo cluster early- 
type galaxies, based on model-independent size measure- 
ments from SDSS imaging data. In both, the relation of size 
and brightness and the CMR, we find noticeable unsteadi- 
nesses between dwarfs and giants. In the former relation 
dwarfs do not fall on the extension of the rather steep se- 
quence of the giants. While the gross trend in the size lumi- 
nosity relation can be explained by light profile shapes be- 
coming steeper for more luminous galaxies, a closer look re- 
veals a statistically significant distinctness in the behaviour 
of faint and bright galaxies. The CMR is continuous over 
the whole range. Yet the observed change in slope and the 
variation of the scatter might hint at more complex reasons 
for this particular behavior than what might be naively ex- 
pected from having one common origin, with the very same 
processes shaping the CMR in the same way. 

But with the comparison to the semianalytic model (not- 
withstanding our rather crude way of selecting early types 
in the SAM) we show that neither of the two findings nec- 
essarily implies a formation by substantially distinct pro- 
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cesses. Instead, the qualitatively similar distributions of the 
model galaxies in the two scaling relations might hint at a 
formation in a cosmological context for the dwarf galaxies, 
hence an origin common to the giant ellipticals. This is in 
accorda nce with previou s claims of no distinction b etween 



them ( Ferrarese et al. 2006t Gavazzi et al. 2005 



Graham & Guzmanll2003l: iMisg eld et al.l l2008h 



It must be mentioned that different approaches can also 
successfully explain the dwarf behavior. For example, ram 
pressu re stripping can re produce the radius brightness re- 



lation (IBoselli et al.l 120081 see also A. Boselli, this issue). 



Given the newly appreciated diversity of dwarf appearances 
we see that as an advantage rather than a shortcoming. 

For future studies, the following procedure seems prom- 
ising: Both relations seem to be r elated in on e or the other 



way with the galaxies' dynamics. iBernardi et al.l d2005) con- 
cluded that the CMR is a result of two more fundamen- 
tal relations: the Faber-Jackson relation and a relation be- 
tween color and velocity dispersion. Given the slope change 
of the Faber-Jackson relation from giants to dwarfs 
(de R iicke et al. ,2005; Matko vic & Guz man, 2005) a change 
of slope of the CMR would actually be expected. And in 
the semianalytic model the sizes are strongly influenced by 
the dynamic feedback. In the model's context this should 
be closely related to the internal dynamics of the galaxies. 
Therefore it is desirable to obtain central velocity disper- 
sions and kinematics for early-type dwarfs in greater num- 
bers. 
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